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ABSTRACT 


Title of thesis: Fine Structure In 3C 120 and 3C 04 

Laura Katherine Hutton* Doctor of Philosophy* 1976 
Tlies Is directed l>y: Professor William C. Erickson 

Seven epochs of very- long-base 1 liie radio Interferometric 
observations of the Seyfert galaxies 3C 120 and 3C 04* at 3*0— cm 
Kavelength using stations at Westford* riassachuse t ts , Goldstone* 
Oaliforiila* Green Bank, West Virginia* and Onsnla* Sweden* have been 
analyzed for.soiirce structure. An algorithm for reconstructing the 
brightness distribution of a spatially confined source from fringe 
amplitude and so-called ‘‘closure phase" data has been deve loped and 
successfuliy applied to ar 1 1 f Ic la 1 ly genera ted tes t data and to data on 
the above-mentioned sources. 

Over the two year time period of observation* 3C 120 ivas observed 
to consist of a double source showing apparent super-re la 1 1 vis t Ic 
expansion and separation velocities. The total flux changes comprising 
one outburst can be attributed to one of these components. 3C 84 showed 
much slower changes * evident ly involving flux density changes in 
indlvldiial stationary components rather than relative motloii. 

Some previously proposed explanations for the behavior of compaat 
radio sources are briefly re-examined In the light of the present maps. 
The possibility that the apparently divergent behavior of the two sources 
studied Is in fact due to variations of a similar mechanism is also 
d iscussed . 
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INTRODUCTION 


"The history of radio Interferometry Is one of steadily Increaslner 
baselines to obtain ever higher resolution" (Bare et ol. 1967)# And 1967 
^ms a good year In the history of radio In ter f erorae try • Early in that 
year, the first interferometers with totally unconnected elements were 
tested, nearly simultaneously In Canada (Broten et ol. 1967a) and In the 
United States (Dare et ol. 1967), and baseline lengths leapt from the 
order of 100 km to virtually the diameter of the Earth. The first 
fringes on a transcontinental baseline (Algonquin Park, Ontario to 
Penticton, British Columbia) were obtained on May 21, 1967 (Cohen 1969), 

and, by July of the same year, angular size measurements on a large 
number of sources were in progress (e. g. Broten et nl. 1967b; Broten et 
al. 1967c; Clark, Cohen, and Jauncey 1967; Clark et al. 1968a; Clark et 
al. 1960b; Kellermann et al. I960; Gubbay e t al* 1969a; Gubbay et al. 
1969b). As the supposed synchrotron self-absorption cutoffs in the radio 
spectra and the rapid time variations had suggested, structure was 
apparent In many compact sources on scales from 1 to 10 milliseconds of 
arc at centimeter wavelengths* 

The early exporiments. Intended to search for Interesting objects, 
provided only very scanty sampling of the fringe visibility function. 

Only general statements could be made about the approximate angular size 
of the brightness distributions. Tlie t heore t lea 1 ly expec ted 
relationships (see e. g. Kellermann and Paullny-Toth I960), between source 
size and radio spectrum were, however, confirmed. 

The first systematic sampling of a source’s visibility function 
came ns a by-product of an attempt In October 1970 to measure 
gravitational light bending by the sun (Ilobertsoii et nl. 19<1). The two 
quasars 5C 273B and 3C 279 were tracked for the entire period of mu tun 1 
visibility for several days, as the latter was occulted by the sun. 
Tell-tale varlntlons in the fringe amplitude ns a fuiietlon of time of day 
revealed that both sources had observable structure (Knight et al. 1971). 
The simplest Interpretation was that of a double source in both cases. 

In February 1971, the experiment i>ms repeated (without the sun), 
revealing changes In the structure of both sources (lOiltney et ol. 1971). 
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Under the snine simplest nssump t Ions , both sources sliowed apparent 
velocities of expansion well above the velocity of light* 

These startling results provided Impetus for the establishment of a 
soiu'ce structure monitoring program on the so-called ‘‘Go Ids tack" (or 
Go Ids tone/Hays tack) interferometer, which began In April 1971 under the 
guise of the "Ounsar Patrol". By May 1973, array e icments were added In 
Vest Virginia and Sweden, Increasing greatly the regions of the (u,v) 
plane which could be sampled. The i>resent study will concern Itself with 
data from these observing sessions. 

In particular, the objects chosen for study are the closest of the 
strong ones, the radio galaxies 3C 04 ( NGC 1275) and 3C 120* In angular 
size, they are among the largest sources that we observed, perhaps for 
that reason, so that a given In* ' m* f erome ter provides th; best relative 
resolution and also the best absolute (as measured in parsecs) resolution 
that It was capable of producing. Fur thei’inore , the distances of the 
radio galaxies are more widely accepted than those of the q[uosars , a 
point which must be considered when computing apparent velocities of 
expansion. 

Chapter II constitutes a description of the experiments themselves 
and the data reduction procedures used. More general mapping techniques 
than the traditional Gauss Ian-component modeling methods proved necessary 
to handle the apparently complex structure o f 3C 04, trspeclally, and to 
make full use of the "closure" method (Rogers et nl. 1974) of extracting 
phase information from VLBI data. An algorithm for reconstructing a 
genei'al coiif Ined source from the data ws developed for this study, and 
is discussed in Chapter III. Chapters IV and V present the results of 
application of these algorithms to the data on 3C 120 and 3C 04. And 
finally. In Chapter VI, conclusions are draim about both the mapping 
algorithm and the sources. 


- 2 - 


Chapter II 


THE QUASAR PATROL 

Tlie experiments that will be considered were conducted between 
April 1972 and October 1973 at a frequency of 7.85 GRz (wavelength 3.8 
cm) using left circular polarization* with some subset of the following 
stations: the 120-ft (36.6-m) antenna at Haystack Observatory ( "H") near 

Westford, Massnchuse t ts i the 2l0-ft (64-m) "Mars" antenna at the 
Goldstone Tracking Station ("G") near Barstow, California) the 140-ft 
(43-ra) telescope at the National Radio Astronomy Observatory ("N") at 
Green Bank, West Virginia) and the 85-ft (25.6 m) antenna at the Onsala 
Space Observatory In Sweden <"S"). An 85-ft (26-m) dish at Fairbanks, 
Alaska ("F") was also running on some of the early sessions, but the 
amplitudes are virtually uncalibrated and are not of much use for the 
analysis of source structure. All data was recorded on the NRAO Mark I 
VLBI system. The sessions ranged from 12-hour Go Ids tack runs to 5-day, 
4-statlon marathons that recorded upwards of 3 tons of magnetic tape. 

This chapter endeavors to describe the observing sessions 
themselves, from plannlsig and scheduling, through set-up and running of 
the stations, and the Initial stages of data reduction and calibration, 
from the point of view of the astronomer Interested In source structure. 
For analysis of some of the same data from the points of view of 
astrometry and geodesy the reader Is referred to Witte Is 1975 and 
Robertson 1975. 

A. THE GEOMETRY OF VERY LONG BASELINES 

Tlie geometry of the Interferometer is of Interest to the astronomer 
for purposes of fringe processing, of course, and because It delermlnos 
the sampling of the visibility function. A summary of the subject Is 
therefore in order. 

The baseline itse 1 f is expressed as a vector which points from 
"atfitlon 1” to “station 2“ » represented as the difference In geocentric 
Cartesian coordinates of the two s ta t Ions expressed In 

llght-ralcroseconds. The fringe spacing at any given moment Is determined 
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by the apparent lengrth of the baseline as seen froni the source (Figure 
II. 1), or, In other vords, projected onto a plane perpendicular to a unit 
vector pointing tovard It. The two-dimensional coordinate ays tern (u,v) 
la defined In this plane, with u Increasing to ihe east and v increasing 
to the norlh, The delay, between the arrival of the same wavefront at 
the two stations, Is Just the projection of the baseline onto the source 
unit vector* Wlien all these quantities are measured In wavelengths, u 
and V represent the two-dimensional spatial frequency that the 
Interferometer samples and the delay Is the total fringe phase. 

As the earth- fixed coordinate system rotates under the sky, the 
projected baseline traces put an ellipse In the (u,v) plane. Thus a 
number of d 1 f f eren t spa t la 1 frequencies may be sampled with the same pair 
of stations. Tlie eccentricity of the ellipse depends on the declination 
of the source, forming a circle for a source at the North Pole and a 
straight line for one on the equator. As a result, the coverage of the 
(u,v) plane Is generally poorer for the more southerly sources. Also 
radio telescopes generally have limit stops, the ultimate limit stop 
being the Earth itself, so that the entire ellipse Is seldom sampled. 

With multiple baselines the coverage Is much improved, but often still Is 
far from Ideal. This sorry situation marks the greatest single obstacle 
to mapping of compact radio sources through VLB I. 

To maximize the coverage. It is desirable to have the various 
baselines in an array as different In length as possible, for minimum 
redundancy. "Crossing points", or places In the Cu,v) plane sampled by 
more than one baseline, can be useful, however. If calibration of each 
baseline Is not possible, crossing points may be used to tie together at 
least the relative scale factors between d 1 f f erent base 1 Ines . 

Before leaving the topic of geometry, consider a triplet of 
stations. The base 1 Ine vectors sum to zero, of course, and so do the 
delays and hence the phases, for a point source. If, however, the radio 
source is resolved to a different degree on the different baselines, Its 
apparent loca t ion may d if f er , and then the phases will not "close". This 
deviation from zero provides another clue to the structure of the source. 

B. BANDWIDTH SYNTHESIS FOR THE ASTRONOMER 

Both geodesy and astrometry depend upon accurate jncasureinents of 
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F Igure 11,1. The pro Jcc ted Hays tack/Onsa In base line 
northerly radio source. As the Far th ro tn tos ♦ 
hnsclliio will trace out on ellipse In the ( ii » v) 
positive to the north and u is pos I live to the 
direction ot Increas Ing r Iglit ascens Ion. 
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the geometric de Iny, Here , however , VLB I is severe ly hampered by its 
bandwidth restrictions, which come about because of the limited data 

capacity of the tape recorders. The sky is full of fringes, so to speak, 

and they all look alike. It Is Impossible to isolate a source in any one 
of them. Fortunately, the necessary wide bandwidth does not have to be 
contiguous; it may be synthesized from an array of narrow pass bands in 
the same way that an aperture may be synthesized from an array of 
antennas . 

In fact, it is instructive from the astronomical point of view to 

consider the bandwidth synthesis system In the (u,v) plane. At each of 

several different frequencies, a single interferometer has a different 
resolution; it becomes. In effect, several Interferometers. The beam 
pattern of the resulting array no longer consists of Identical fringes. 
The " whi te— 1 Ight fringe" has been, Isolated to an extent which depends on 
the maximum frequency spacing, and with an ambiguity which depends on the 
minimum spacing. The geodesist usually prefers to think of the beam 
pattern in the delay domain and call it a "delay resolution function". 

He measures the geometric delay as the rate of change of fringe phase 
with frequency, rather than with resolution. Tla thema t lea 1 ly, however, he 
does precisely what nnastronoraerdoesto locate a point source with a 
synthesis array. (For the gory details of theory and practice see Rogers 
1970, Hlnteregger 1972 and Vhltney 1974.) 

Thus the typical Quasar Patrol experiment Is spent switching in 
rapid success Ion ( 0 . 2 second intervals) be tween 3 or 6 narrow passbands , 
each 360 kHz wide, spaced over about 50 MHz near the "center frequegtey" 
at 7850 MHz. The fringes are processed in each "channel" separately and, 
if poss ib le , averaged coherently with the phase reduced to this center 
frequency. Ideally, there is no loss o f sens 1 t ivl ty over a non-swltched 
observation. Ins trumen ta I d ispers ion must be calibrated by Injecting a 
coheren t s igna I a t the beginning o f the rece i ver cha in a t each s ta t Ion 
and processing these signals as If they were fringes. The "phase 
calibrator" signal is injected cont liiuous ly, but is weak enough not to 
affect the system temperature, and, since it correlates at zero fringe 
rate. It does not ordinarily Interfere with the detection of the true 
fringes. In the event of malfunction of the phase ca 1 ibra tor , the 
fringes of a strong radio source can be used Instead, assuming that the 
Ins trumen ta 1 d Ispers Ion is the same for the o their observe t Ions as we 1 1 . 

If the Instrumental d Ispers ion Inf orraa tlon Is lost, the amplitudes from 
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the individual channels must be averagfed Incoherent ly^ with a 
corresponding; loss in sensitivity. 

It may one day be possible to use the extra ( v) lane coverage 
gained by frequency switching for the analysis of source structure. With 
the presen t s ys te in » the f r Inge amp 1 1 1 udes and phases do no t change 
appreciably over the less “ than"”one "•percent spanned baud widths and the 
increased sensitivity gained from the coherent averaging Is consideriid to 
be more important. Doth greater recording and spanned bandwldths arc 
planned, however, for the future Mark 111 system. It may then be 
advantageous to treat each frequency channel separately in source 
structure analysis. 

C. SCHEDULING 

From the source structure point of view, (u,v) ““plane coverage is 
the main concern; for any one source, it is des irabe to make measurements 
over as wide a spread in hour angle as possible. Given unlimited 
telescope time, or interest in only a few sources, and no geodetic 
requirements, the ideal schedule would have all the stations following 
one source at a time through its entire observable hour angle range. 

Only then is there a chance of obtaining the maximum amount of 
Information about a visibility function that a given VLBl array can 
provide. This sort o f schedu 1 Ing, ho wever , was tes observing time since 
some telescopes have to spend time waiting for sources to come up 
some whe r e e 1 s e . 

Fortunately, most of the sources of concern to this program do not 
require such drastic measures. The visibility function can be expected 
to vary smoothly for a period of time dependent on the angnlar s ize o f 
the brightness distribution. W 1 th the except Ion o f 3C 84 , the sources 
are small enough to allow gaps in the data on the order of an hour. llie 
practical solution to limited observing time, then. Is to move the 
antennas back and forth between a few sources In the sky, returning to 
each often enough so that the visibility function Is "reasonably well 
determined". It should be noted that this sort of scheduling Isnlso 
wasteful, In that large amounts o f t 1 me are spent moving the telescopes. 

I t does , however, better correspond to the geodesist’s scheduling 
precept! the antennas should be eye ling between at leas t three sources at 
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once . 


A fi^eneral kno>v^ledge of the expected visibility function for each 
source Is necessary to make a schedule « as well as knowledge of the time 
scale of Its variability (Vhich sources must be included in every 
exper Iraent?) . In addition to source structure, astrometry, and geodesy 
re qiulretnen ts , there are other headaches for the scheduler* For example, 
differential experiments (see Witte Is 1975) tend to monopolize large 
blocks of time, causing gaps In the data for other sources* And there is 
usually a list of "trial" sources to be allotted a few runs each during 
the session. Also, it takes any normal human 30 seconds to change a 
tape, not counting rewind time. Some te Icscopes , typically Goldstone, 
are not available for the whole experiment. Different telescopes move at 
different rates about different axes. (A particularly tricky situation 
occurs when two sources near the same right ascension transit at 
Goldstone, one on the north side of the zenith and one on the south. 

This problem occurs with 3C 84 and 3G 1,20.) Telescope limit s tops are 
ail different. Cables must not be wrapped around the azimuth axis by 
more than their length. Extra time must be allowed for "peaking up" and 
for radlometry. Previous to the May 1974 session, Onsala had to be 
allowed two hours at regular Intervals to replenish the cryogenic system. 
And so on, ad nauseum. 

For these reasons, and more, efficiently scheduling the Quasar 
Patrol array has become an art which has never been completely 
computerized. The first appi’oxima t Ion to a schedule, In fact, is often 
constructed with a celestial globe in hand, along with tables of rise and 

set times. Blocks of schedule may be lifted from previous experiments. 

One day’s schedule should, as far as possible, be complementary to the 
previous day’s In sidereal time, to fill in the gaps In the sampling. 
Extra runs may be Inserted for the f as t~"ino vl ng stations while the slow 
ones are still changing source. The proposed schedule is then checked, 
to make sure that all the stations can keep up with it, by an ingenious 
nol-so- 1 1 1 tie computei* program (Vandenberg 1974). After a suitable 
number of iterations, the scheduler may be blessed with a "final" 
schedule, which Is sure to be changed at the last minute over the 

telephone. Tlie same computer program then prints multiple copies of 

f 11 1- In- the-b lank observing logs for each station. 
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D. SETl'ING UP AND RUNNING THE EXPERIMENT 


As the appointed time approaches, the experimenters disperse 
themselves atnong^ the various stations and begin setting up for the 
experiment. Since the Mark I VLBl regime represents a major perturbation 
on the normal operation at most stations, this process can take as much 
as a day or more. 

The hardware re qulreraents for both astronomical and geodetic VLBl 
have been described so many times (see e. g. Hlnteregger 1972, Whitney 
1974) that It hardly seems necessary to add more words. A block diagram 
of the Green Bank station, however, Is Included (Figure 1 1. 2). The other 
stations, and even Green Bank on occasion, may have differed in detail 
from this example, but the general principles must remain the same or 
there are no fringes. 

A few restrictions were put on the source structure effort by the 
requirements of the others, and benefits were reaped as well. One of the 
most far reaching in effect was the choice of the NRAO Mark I recording 
system, over the Mark II system. The Mark I’s narrower recording 
bandwidth (360 kHz) restricted the sensitivity and limited study to the 
brightest dozen or so compact radio sources. Tliree advantages, however, 
emerged from its use. The first was that a great deal of care had been 
taken to account for all losses In correlation, so that it was possible 
to calibrate the fringe amplitude scale without appealing to unresolved 
sources, which may or may not exist. The second wind fa 1 1 came from the 
care the geodes 1 s Is have taken to preserve the meaning of fringe phase, 
making It re la t l ve ly easy to construct the so-called closure phase, on 
which much of the fo 1 lowing aria lys Is depends. Lastly, the Mark 1 system 

Is more reliable than the Mark M. 

Of particular Interest here are the aspects of the hardware which 
affected t he f r Inge amp 1 I tude ca I Ibra t Ion. Discussion of these , ho wever , 
is best understood In the context of the data process Ing and wi I V be 
postponed to a section on that topic, 

The actual rimiilng of a VLBl session can bo, with Inck, extreincly 

dull. When everything la working-, the dut lea o f the observer cons Is t of 
changing the tapes, cleaning the heads , making the tape labels, keeping 
the logs , making sure the antenna opera tor po Ints the antenna at the 
right place at the right time and takes the calibration data, and other 
even loss exciting chores. The long periods of boredom are some times 
























broken, however, by momeula of intense panic, diirlngr which data Is 
usually beliiff lost. It then bocornes the observer’s responsibility to 
find and fix the problem, or find someone who can. VLBI suffers from the 
added complication that, although wany checks on part or most of the 
system are run routinely, no overall check is available in real time, if 
fringes are not later found, 1 t is not always possible to reconstruct 
what happened. Constant meticulous attention to detail and large amounts 
of faith are required of a VI-.BI experimenter, is a tribute to the 

skill and dedication of the (large) Quasar Patrol team that so little 
data has actually been lost# 

At the end of each session, the enormous pile of tape boxes outside 
the door at each station was shipped off to Haystack for processing. 


E. DATA REDUCTION 

In the case of VLBI, the terra “data reduction" inay be taken quite 
literally. In the first of two steps, the program “VLBH" compresses on 
the order of 100 runs onto one output tape, which includes the complex 
correlation coefficients of each 0.2 second "record'* of data for seven 
different delays. Tlie program "VLBI2” fur tlier reduces the tape tonnage 
to 500 runs per output tape, by searching for and recording only the 
maximum (complex) correlation for each record. Tlie 3 tons typically boll 
down to 4 or 5 VLBI2 output tapes. The Inner workings of both these 
algorithms have been descr Ibed in great de ta i 1 by Ilin teregger ( 1972) and 
Whitney ( 1974) In their theses, and In summary form by Vandenberg ( 1974) , 
and will not be discussed further here. However, some description of the 
output format and of the correc t ions made to the fringe amplitudes Is 
relevant to the calibration problem(s). 

At this point, the data consist of, for each record, the complex 
correlation coefficient at the fringe ra te and delay for which the fringe 
amplitude was a maxlinum. Corrections have been made for the effects of 
Infinite clipping, nonzero residual fringe rate, etc., so that the 
amplitude gives the fraction of the recorded noise which was correlated. 
The Mark I system and the data processing programs were checked (see 
Wlttels 1975) by simultaneously recording tapes Including independent 
no ise wi th a known amount of corre la ted no Ise added , Tlie result Ing 
correlations were alivays within 6 percent of that expected over a large 
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range of fr inge ra te and aiti|> 1 1 1 ude > 

For an Ideal, s Ingle*- f re <luency run, the complex amplitudes could be 
simply averaged over the whole lengtit of the tapes to Increase the 
Integration time* Because of the bandwidth syntheslsi howeve^, the 
records are not alt taken at the same frequency. Instrumental dispersion 
and the change In resolution with frequency cause phase rotation. The 
practice, In the mu 1 1 1 f re quency case, Is to rotate each channel to the 
center frequency, using the phase ca 1 Ibra tor to determine the 
Instrumental dispersion and the observed group delay to correct for the 
resolution difference. It should then be possible to average the complex 
correlation coefficients coherently through the whole tape. Problems, 
however, arose In practice. Some caused only headaches and others 
significant loss of data. 

Because we are forced to work with the amplitude of the fringes, 
irrespective of phase, noise causes a positive correlation even when 
there are no fringes. Wliltney ( 1974) has worked out a theoretical 
estimate of the expected correlation In the no-slgnal case (for clipped, 
band-limited data sampled at the Nyqulst rate): 

1 jC 

' A 

where Is llie bandwidth and >&t the integration time. Units are, again, 
fraction of the recorded noise. This figure serves as an estimate of the 
uncertainty In the es 1 1 nia t ion o f the correlation coefficient In the case 
where there are fringes* With 360 kHz bandwidth and 3 minute coherent 
avejrngcs, the conventional 5 sigma detection limit becomes about 0.0007, 
or approximately 0.3 Jy on baselines Involving the Go Ids tone antenna* 

The correspond Ing figure on HN is about 0.5 Jy, and for the Sweden 
baselines higher still. 

Where there are fringes, it Is possible to allow statistically for 
the add 1 1 Iona 1 corre la t Ion that the noise causes: 

z 1 ' ^ 

e - p, - 0? 

where p Is a better estimate of the correlation than p^, which VLBI2 
estimates. This correction Is made to all the data, but Is only 
significant for weak fringes, with correlation coefficients less that 
0.0020* When there are no fringes, even the corrected fringe amplitudes 
are still positive, although small, and can systematically bias some of 
the amplitudes upward and perhaps affect the source maps. 
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If the system temperature is not the same in all the frequency 
channe is T then corrections wouid have to be made to the amplitudes as 
well as the phases. Ve have founds however, that variations are 
negligible except when one channel la very near the edge of the receiver 
passband at some station. Rather than try to correct the amplitudes, we 
threw the data from the bad channels away, 

Another, more recurrent problem, was phase noise due to somewhat 
less than coherent local oscillators. If a coherent average is mode for 
a time longer than the time scale of the phase variations, the amplitude 
is systematically reduced by an unknown factor. Shorter coherent 
averages (normally 100 records or 20 seconds) were then used, with an 
ensuing sacrifice In Integration time. The amplitudes from the 
individual 100 record "segments'* were averaged Incoherently to produce a 
whole tape average with a larger error bar. 

Other equipment problems caused trouble, as well. Beginning with 
the August 1973 session a recurrent tape recorder problem at Green Bank 
and Sweden caused Incredible numbers of parity errors and significant and 
variable loss of correlation. The error was almost completely 
compensated for In processing beginning In March 1974, but the two 
week— long, 4— station August and October 1973 sessions were left In 
shambles, most of the data tapes having been re-used before the cure was 
discovered. Large amounts of bad and doubtful amplitude data wore thrown 
away, and many more were re— averaged using only records without parity 
errors. As a result, much of the data from these experiments Is 
considerably noisier than the other sessions. The August session had 
calibration trouble as well, as we shall see shortly, so the loss was no t 
as bad as It could have been. 

During the March 1974 experiment an electronic problem at Green 
Bank caused a more or less uniform loss of correlation, the amount of 
which had to be determined empirically. This correction depends on the 
system temperature, but, to a close approximation, It can be assumed to 
be constant for any one source. 

At this point, we again summarize the state of the amplitude data, 
now considerably "reduced". With luck, we have an estimate of the 
fraction of the recorded noise which was correlated, along with an error 
bar. These data must now be merged with the station radlometry to 
translate the amplitudes Into units of flux density. 
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OElGKAIiPA®® 
OP POOR QUALm 


F. EADIOflETllV 


To norma llzo the frlng^en. It Is necessary to know the fractloiial 
contribution of the radio source to the total system noise temperature at 
each station. In practice, both the system temperature (T^) and the 
contribution to It due to the source ( T^) are either measured or 
calculated individually. Both values can «;enerally be expected to vary 
Kith antenna elevation, and with the weather, so that a more or less 
continuous record is virtually mandatory for a g^ood calibration. 
Unfortunately, such a cont Inuous record has not always been possible and 
it has been necessary to Improvise some of the ca 1 Ibra t ion. 

The simplest radlometry scheme was that employed at both Coldstone 
and Haystack, where the (broadband) sys tern temperature was measured using 
a nolse-addlng radiometer system. Measurements both on and off source 
were repeated, if possible, before and/or after each VLBI tape. This 
system Is par t Icular ly e f f ec t i ve at Coldstone, where the system 
temperature is very low, the effective aperture is very large, and the 
weather Is usually very good. The Coldstone radlometry data has also 
been used whenever possible for determining total source fluxes for the 
times of the experiments. Tlie source fluxes were scaled to total power 
observations of DR 2 1 and 3C 274 (Dent 1971, Baars and Hartsuljker 1972). 
The Hays tack antenna is smaller, however, and the receiver and atmosphere 
noisier. It was not always possible to “see " the weaker sources in total 
power. The usual remedy, in such cases, was to compute the source 
temperatures based on an eff Ic leucy-versus-e leva t Ion curve determined 
from a stronger source and the source flux as measured at Coldstone. 

At Green Bank, the sys tern tempera ture was monitored by a manual 
version of the no Ise-add ing rad lome ter scheme, with a measurement taken 
for every tape. The source temperatures were measured, however, in a 
beam-switching mode, whl ch e 1 Imlna ted the effects o f , among o ther things, 
the beautiful West Virginia clouds. As a result, the Green Bank 
calibration was second only to Go Ids tone ’ s and was sometimes used for 
to ta 1 fluxes . 

As far as calibration is concerned, the "problem child" was 
definitely Onsala. The antenna Itself Is relatively small and wns not 
Intended for use at 7.8 GHz; the aperture efficiency was bare ly 15 
percent. On top of this, no successful effort a t cal Ibra t Ion was made 
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until near the end of the experiment series. One attempt did occur In 
August 1973, when several monsurements were made of Cassiopeia A to 
determine the antenna efficiency. All the data were recorded on a strip 
cliart, which unfortunately had to be cut Into a number of sections to get 
It off the malfunctioning recorder. The result looked loo much like 
trash and was Inadvertently thro>\m away before It could be analyzed. 

System temperatures for each observation, with the exception of August 
1973, were measured with a “Y-factor“ system, wherv^ a variable attenuator 
Is used to compensate for an Increase In system temperature due to a 
noise source* The attenuator setting (the ’'Y-factor“) was recorded In 
the log. For sessions prior to May 1974, all the source temperatures 
were computed based on an efficiency of 12 percent as measured In May 
1973, on the assumption that pointing was perfect and that no elevation 
dependence was present. Tlie former Is probably not a bad assumption, 
except In high winds, whereas the second could cause large systematic 
errors. In May 1974, a new cryogenic system was Installed, and the noise 
source temperature was changed and never measured accurately. A 
concerted effort was made to calibrate the new system In July 1974; much 
to everyone^ s surprise, 3C 84 (at 60 Jy) was dlscernable In the records. 
(The Y“factor attenuator was used to compensate for the Increase tn the 
no Ise tempera ture due to 3G 84, and the reading recorded.) During this 
particular session, the system temperature varied from about 105 K. near 
zenith to 120 K. a t 15 degrees elevation, and Increased rapidly below 
that, Tlie scatter on the eflciency measurements was close to 20 percent; 
within that range, there was no evidence for any variation with 
elevation. The average value of 0.03 K/Jy, or about 15 percent, was 
adopted for both May and July 1974, even though some Inconsistencies 
result (Wittels 1975). The scatter Is much less in the correlation 
coefficients, Indicating that most of the noise originated in the 
radloraetry methods, rather than, say, pointing. The upshot of this 
discussion is Hint rather large systeinallc errors probably exist In the 
Onsola calibration. The overall scale could easily bo off by 20 percent. 
In addition, a "skew" In the calibration scale, as a function of Onsa la 
hour angle, might exist due to neglect of the elevation effects 
(especially In August 1973, when the system temperature had to be assumed 
constant). An intuitive estimate of the larges t poss lb le systematic 
error leads to a figure In the 20 to 30 percent range. 

Estimates of the error bars for the rad tome try have generally not 
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been too rigorous. In cnscs where tho temperatures were inonsured or the 
antenna parameters well determined (In other words, II, N, and G) , they 
were estimated from the scatter, which seems to bo largely duo to 
atmospheric effects, and Is therefore time variable, Wlien fringes are 
strong, most of the error In tho final normalized fringe amplitude comes 
from the radlomotry, so It Is a stroke of luck that our best calibrated 
antennas determine the shortest baselines. Radlometry errors for Onsala 
were estimated from uncertainties In the parameters used to calculate the 
temperatures. Fringes wore seldom strong on Onsala baselines and most of 
the error In the normalized amplitudes came from the uncertainty In the 
correlation coefficients. 

Radlometry data aiid error bars having been determined by one means 
or another, It remained only to normalize the correlation coefficients to 
units of correlated flux density 


Sc’ SpJ 


Ti.1 Toa 


Here S Is the total flux and the fringe amplitude in Janskys . It has 
been customary to work with the fringe amplitude in these units because 
the maps may then be computed in brightness temperature units If desired 
Error bars for the fringe amplitudes were computed from the correlation 
coefficient and rad lorae try error estimates, assuming all the 
contributions to be Independent. In most cases the computed error bars 
agree reasonably well with the amount of scatter in the data, a fact 
which inspires confidence In the estimation procedure. 


G. CLOSURE PHASE 

To measure a source br ightness dls tr ibut ion uniquely, we must know 
the complex visibility function. Amplitude data a lone , no ma t ter how 
much of the (u,v) plane it samples, cannot produce a unique so lut ion. 

VLBl Is not generally able to measure the fringe phase because the 
relative phase of the local oscillators at the two stations as a function 
of time, and atmospheric phase shifts as a function of time, r”e unknown. 
We are therefore forced to "make do" with the next best thing, the 
so-called "closure phase". Tlie Idea was first used by Jenn iso n C 1950) 
for early radio link Interferometers, >vhlch were no t phase s tab le , and 

has been applied in Its present form to VTjBI data by Rogers et al* 
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(1974). The closure phase Is simply the sum of the observed fringe 
phase, Including all of the error terms, summed around the triangle of 
baselines formed by three stations. All phase error common to only one 
station, including local oscillator and atmospheric phase errors, add 
twice with opposite signs and cancel. Ve are left with the sum of the 
visibility phases on the three baselines. 

One minor complication Is Introduced because a triplet of stations 
(e. g. "HNS") Is gencirally processed in the order HN, HS, NS* "Baseline 
2" (HS) is backwards and the sign of Its phase must be reversed for the 
sum. 

Another minor complication comes from the fact that the earth 
rotates between the times of arrival of a wavefront at the different 
stations. Hence, by normal processing conventions, the delay (and lienee 
the phase) on baseline 3 has been computed for a slightly different a 
priori time than on the other baselines. The baseline 3 phase must be 
rotated back In time to correspond to the arrival time at station 1 ( in 
this case H) . 

Frequency switching remains a minor comp I ico t Ion as long as the 
phase calibrator Is working well* In that case, the Individual channel 
phases cun be rotated back to the center frequency, averaged, and summed 
to obtain a closure pliase appropriate to that frequency. If the phase 
calibration Is poor, this procedure will Introduce extra noise Into the 
closure numbers but will not produce any systematic effect. In no case 
is any of our data affected by more than about 20 degrees of phase due to 
poor phase calibration. 

Phase noise does not affect closure phase, so long as the 
Individual baseline coherent averaging is done over a time period shorter 
than the time scale of variation. If the coherent averages are longer, 
however, noise can be introduced, but again there are no systematic 
effects. To avoid the problem, closure phases were constructed for each 
100-record segment and these were averaged Incolierent ly for the whole 
tape . 

Closure plinse for runs severely affected by parity errors had to be 
recomputed, deleting all the bad records. 

The error bars which were computed for the closure phase data 
points include only the uncertainty ((T<^ ) in the phase estimation Itself 
(see Whitney 1974) 5 fj^ j 
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The closure phase error bar Is taken to be a Quadratic sum of three of 
these, one for each baseline. This is appropriate when the fr Ingres are 
strong* Tlie size of the error bars Increases rapidly with decreasing 
correlation, meaning that the closure phase is noisy when fringes are 
weak on any one of the three baselines. The additional noise due to poor 
phase calibration and phase noise has not been included; in some oases, 
then, the size of the error bars may be too low. Again, however, the 
error es t ima tes • agree reasonably well with the scatter in the data. 

It should be pointed out that for a 4-statlon run, only three of 
the closure phase data points are Independent. The fourth is a linear 
combination of the other three, and adds no new information to the 
so lut ion. 

H. SUTIMARY OF EXPERIMENTS 

For convenience, a table (Table II. 1) 1 is t ing a 1 1 the Quasar Patrol 

experiments has been prepared. Sessions have been Included in the table 
whether or not they contain data on the sources under consideration here. 
Tlie station abbreviations (H = Haystack, G = Go Ids tone, N = NRAO, S = 
Sweden, F = Fairbanks, Alaska) will be used throughout this thesis. Any 
problems which could have affected the accuracy of the data, such as 
problems with the weather, the frequency standards, or any of the other 
equipment, are listed. Unless otherwise stated, hydrogen maser frequency 
standards were used at all the observa toi^ ies . 
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Table II. 1 A List 

of Ctunsar 

Pa tro 1 

Experiments. 

Da tes 

Slat Ions 

Problems 

14, 15 April 1972 

II, G 



9, 10 JIhy 1972 

H, G, 

F 

Rain at H 

29, 30 fUiy 1972 

H. G, 

F 


3 June 1972 

II, F 



6, 7 June 1972 

II, G, 

F 


27, 20 June 1972 

II, G 




3, 4 July 1972 

II, G 

Not frequency switched 

29, 30 August 1972 

n, G 


7 November 1972 

H, G 


4, 5 February 1973 

II, G 


30, 31 March 1973 

II, G 


17-19 May 1973 

II, G, S 

Rb standard at G 

22, 23 riay 1973 

II, s 


10-14 August 1973 
Cs standard at G 
No S radioinetry 

II, N, G, S 

Parity errors at N and S 

12-16 October 1973 
Rb standard at G 

II, N, G, S 

Parity errors at N 

22, 23 January 1974 

II, G, S 

Scheduling problems 

4-0 March 1974 

II, N, G, S 

Equipment problems at N 

29 Apr 11-3 May 1974 

II, W, G, S 

Rain at H 

0-12 July 1974 

II, N, G, S 


27 , 20 Oc tober 1974 

II, N, G 


15-17 January 1975 

II, N, G 


20 May 1975 

II, G 


24 Angus t 1975 

II, G 
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Chapter III 


DATA ANALYSIS 


Given seven boxes of cards, containing the amplitudes and closure 
phases and error estimates, the observer’s next chore Is to translate 
these Into source maps. This Is precisely the problem over which 
VLBl’ers have torn out hair ever since there has been enough data to even 
consider the possibility. It Is not on easy one, for basically two 
reasons: lack of phase Information and poor sampling In the (u.v)-plane. 

Without any form of phase data, the solution for the brightness 
distribution Is ambiguous. There are at least two solutions, namely the 
correct one and Its mirror linage, and probably many more. The use of 
closure phase provides at least a partial solution to the problem. Even 
with It, however, VLBI does not become true aperture synthesis. Because 
the phases on the three baselines must somehow be separated. It Is not 
possible to work directly with the complex visibility function? the 
algorithms roust be more devious. 

Closure phase makes possible a better description of the visibility 
function at the points sampled, but It does nothing to Improve the 
sampling. Only more and/or bettor distributed baselines can do that. 
Between the case where sampling Is complete and the one where there Is 
clearly not enough data, there lies the hazy region of "almost enough. " 
VLBI with four or five stations, which Is common enough now but 
constitutes a tremendous effort, lies In this hazy realm. Generally 
speaking, It Is necessary to exchange some a priori Information about the 
source for the missing data. A common form for this information is 
assumptions about the number and shape of the Individual source 
components (see e. g. Purcell 1974). It mayor may no t be poss lb le to 

represent a given set of data In this fashion, since sources seldom have 

true Gaussian brightness distributions. In this analysis, the 
assumptions have been made instead about the angular s Ize of the source, 
the goal being to find the smallest source which can represent the data 
rather than the one with the least components. Before the mapping 
algorithm Is described In detail, however, a short digress Ion on the 
interpretation of closure phase is in order. 
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A. ON THE VISUALIZATION OF CLOSUHE PHASE 


It Is eosy to see that the closure phase for bji unresolved source 
Is identically zero, irrespective of the source ’ s locn t ion in the sky. 
Tiierefore, closure phase is of no direct use in astrometry. 

A symmetric source of any typo produces a similar effect, since 
such a brightness distribution has a real Fourier transform. If all 
three visibility phases are equal to zero, the closure phase mus t be 
also. The real visibllty function, however, may be negative, which we 
Interpret as a phase of plus or minus 180 degrees. Then the closure 
phase will be 0 degrees when an even number (including zero) of the 
baselines show a phase of i80 degrees, and will be 180 degrees when an 
odd number do. As a function of time, one would expect the closure phase 
on a given triplet to a 1 terna te abrupt ly be tween 0 degrees and 180 
degrees as the amplitudes on the various baselines pass through sharp 
nul Is . 


Any deviation from 0 degrees or 180 degrees indicates that the 
brightness distribution Is not symrae tr Ic , and also indicates the 
direction of the asymmetry. Tills becomes clear when we recall tliat 
inverting the brightness d is tribut ion reverses the sign of the vis lb i 1 1 ty 
phases, but leaves the amplitude unaffected. The closure phase therefore 

determines the orientation of the source, a p lece of Informa t ion which is 
not available from amplitude data alone. 

As a simple example, let us consider an equal point double 
brightness d is tr Ibut ion on the equator, to be observed by throe east/west 
baselines. Its visibility f unc t ion wou Id look 1 ike F i gure Ill.l.a. In 
the region where the amplitude is negative, we measure a positive 
amplitude with a phase of iOO degrees. At t Ime A, the beginning of the 
track, all baselines are within the first null and have 0 degree 
Visibility phase; hence , the closure phase is 0 degrees also. At time B, 
however, baseline 1 (the longest) has crossed the null and its phase is 
now 180 degrees, while the other two phases are still 0 degrees. At time 
C, ail the phases have returned to ©degrees. The behavior of the 
closure phase as a function of time would look 1 ike Figure 111.1. b. Hie 
same "exciirs ion from zero" effect can also be produced by a tr iplet 
containing two long baselines which cross the same null one after the 
other. In the case of equal point doubles, the closure phases arc 
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relatively easy to visualize. If the source components were not points* 
but had the same shape, the source remains symmetric and the phase 
behavlorlssimilar. 

If the po in t components were une qua 1 in flux, the sharp nul Is wou Id 
disappear and the phase Jumps would be more gradual and would not reach 
180 degrees . If the same array of eas t/wes t In ter for o me ters were to 
observe such a source, the closure phase jumps would also be more gradual 
(see Figure III.l.c). It is easy to see how quickly the closure phase 
curves become complicated when all the baselines are crossing minima? 

It is Instructive to think of a double source with extended 
components as if It were a point double with each component having a 
power equal to the correlated flux the true component would have at the 
current reso lu t ion. We have, in e f f ec t , a po int doub le whose f lux 
densities change as a function of resolution. This Is strictly true if 
the visibility function for each component has 0 degree phase, for 
example, If it is symmetric and is being observed inside the first null* 
Otherwise, the component positions as well as their powers are functions 
of reso lut Ion . 

Although this system of graph-paper Fourier analysis may seem a bit 
artificial, the author has found it very useful in deciding what to 
expect in analysis of a given source. It has saved a substantial amount 
of computer time, if nothing else. 

B. NAPPING 

The more or less traditional method of source s true ture ana lys Is in 
VLBI is model-fitting (see e. g. Purcell 1974). In this technique , the 
brightness distribution Is assumed to consist of a small number of 
individual circular or e 1 1 Ipt ical Gauss Ian components. The parameters of 
the components, relative positions, fluxes and sizes, etc., are then fit 
to amplitude data by some sort of Iterative leas t-s quares technl que . The 
same process can be app lied with c los ure phase da t a , used s Imul t aneo us 1 y 
with the amplitude data, with some Improvement. The closure phase 
de t ermines the or ten tat ion o f the brightness d is tr I but Ion, and s Imply 
adds more data. Mode I- f 1 1 1 Ing, e ven wl th closure phase, however, has its 
difficulties. For one thing, the Ini tial guess at the model from which 
the Iterations begin must be quite c lose , a fact which makes the 


- 23 - 



procedure frustrating^, at best* In addition, the models are confined to 
a relatively small number of simply sliaped components. The actual source 
may or may not be representable to within the errors by the chosen model. 
Furthermore, except for very simple cases. It is virtually impossible to 
evaluate the deg^ree of uniqueness of a model, once it Is found. 

Ve would like to hove a method of producing a map or maps of the 
brigrhtness distribution. A map, in this context, would refer to any 
regular grid of hypothetical point sources whose flux densities are such 
that the observed visibility amplitudes and (closure) phases are 
reproduced. Let us assume, for the moment, that the source has been 
completely resolved on the longest baseline ( ; the fringe amplitude 
is zero there and beyond. Then the well-known sampling theorem of 
Fourier transforms (see Bracewe 1 i 1963) tells us that such a grid 
representation is possible and that the spacing of the grid should be 

Larger grid spac Ings will not reproduce the data and smaller 
spacings extract no new Information. The map grid may be as large as 
desired, but It need only be as large as the brightness distribution. 

The technique adopted here has been to least squares fit the flux 
densities of these hypothetical point sources to all the data, amplitude 
and closure phase, simultaneously. A little thought will show that this 
is equivalent to fitting for the coefficients of the (complex) Fourier 
series coefficients of the visibility function. Unfortunately, because 
the amplitude and closure phase are nonlinear functions of these 
coefficients, Iterative methods are necessary. Specifically, a gradient 
expansion search for minima in space , ''Mar quard t ' s algorithm", has 
been employed. Tlie solution can be expected to be we 1 1 de ter mined 
provided that the brightness distribution Is restricted to be small 
enough that most of the correspond Ing sampling cells In the (u,v) plane 
contain data points* If this requirement Is no t sa t Is f led , the 
(u,v)-plane coverage is deemed insufficient to map a source of the 
required size. 

Comp 1 lea t ions ar ise in cases where the visibility f unc t Ion is no t 
yet zero on the longest baseline (see Figure III. 2). This is the case 
for vlr t ua 1 1 y a 1 1 Uuasar Patrol observations since the sources are picked 
to show fringes on all the baselines. Some sort of spatial frequency 
filter must be applied to force the visibility function to zero beyond 

or in other words , the source brightness distribution must first be 
smoothed to a resolution appropriate for the antenna array before being 
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The sampllnff of a one-d Imens loiial source whose 
fuiictloit docs not go to zero beyond the innxliimin observed 
Power is nllnsed Into the observed part of the visibility 
row the adjacent Image. 
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F Igure I I 1 . 2 . b . Ono-d line ns Iona 1 convolving function and spaclal 
frcfluency filter. 
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Figure III. 2 .C. Smoo thed brightness d Is tr Ibut Ion and filtered visibility 
f unc 1 1 on. 
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sampled. The sinoo thins function Is arbitrary, but If the necessary srld 
Is to bo kept as small as possible there can be no side lobes. The 
function that has been chosen Is the "rectangular tr langle func t Ion" 
whose Fourier transform Is 

This choice Is made over an elliptical function on the basis of computer 
time considerations; the Fourier transform of the "cone function" 

Involves Integ^rals of Bessel functions. The observed fringe amplitudes 
are multiplied by T(u,v) before fitting. We are guaranteed that at most 
one more grid row will be necessary on each s Ide of the brightness 
distribution, because of the convolving function. The maximum error 
possible due to overlap with the next image In the (u,v) plane is about 5 
percent and if the source is well resolved it should be much less. In 
order to avoid setting any hard-won data points entirely to '.pero, and 
hence jeopardizing the least-squares solution by doing the same to the 
error bars, has been set 10 percent beyond the highest resolution 

actually sampled. 

In exchange for such a "clean beam", one sacrifices resolution. In 
this case, however, where we have very little knowledge of the source 
structure to begin with, the absence of confusing side lobes should be 
cons ide red eno ugh of a b less ing to Jus tlfy the loss in re solution. 

Fitting in the sky domain in this way allows the enforcement of 
another physical restriction: the source br Ightness must be everywhere 

positive. The volume of parameter space which must be searched is 
thereby reduced by a factor of 2 , where N is the number of grid points. 

Practice has shown that this procedure Is usually necessary. Evidently 
the noisy and Incomp le te da ta d Is tr Ibut Ion can a I low unphys lea 1 solutions 
to have lower values of than the desired one. The convergence of the 

algorithm Is slowed in the process. 

The use of any sort of i tera 1 1 ve teclin 1 que always brings up the 
question of initial parame ter vn lues . Here the gr Idd ing ine t hod shows its 
character as the model fitting techn 1 que that it is. There is no a 
priori way of knowing what initial grid will a llow convergence along a 
"clear path" in space to the deepest minimum. T1 rough experience , 

however , It has been found that a de ta I led in 1 1 ia 1 approxlma t Ion Is 
usually not necessary. It has generally been adequate to begin with all 
the parameters zero except for one point, which Is set to some small 
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value. If the point haa been selected so that the source will fit on the 
grid with a peak in the brlglitncss at that position, the solution will 
generally converge, A number of trials may be necessary to find such a 
starting point, but the approach Is systematic and requires no detailed 
Knowledge of the brightness distribution other than Its overall size. It 
Is also possible, however, to Input a better Initial approximation. For 
sessions with relatively little data on a given source. It Is sometimes 
convenient to begin the Iterations with a sonrce map from some other, 
better determined, set of data (or from any other suitable form of lucky 
g^iiess Ingr) . 

For presentation purposes, the sky grids wer^ all convolved with 
the transform of T(u,v) and formed Into contour plots. 

C. ERROR ANALYSIS 

We now assume, for the moment, that the form of the fitting 
function alloivs a reproduction of the data for some parameter values. In 
other words, the data contain no systematic errors and the grid la large 
enough. Three considerations then arise In the error analysis. Assuming 
that a minimum has been found, we need to evaluate the purely statistical 
effect on the solution of the error bars on the data points, and also the 
intercorre la t ion between the parameters. We also need to know something 
about uniqueness; how many other similar minima exist within reasonable 
value ranges of the parameters? 

These different effects will be evaluated as we 1 1 as possible in 
this section, but first a fairly detailed description of the 
least-squares f 1 t t ing technique la necessary. The algorithm works In 
parameter space, which has as many d linens ions as there are fitting 
parameters. The function or the mean square deviation of the 

computed curve from the data, Is a continuous function of position in 
parameter space. A typical minimum should be somewhat e 1 1 ipso Ida I in 

shape near the "bo t tom" . Inside the minimum, may be evaluated at 

various points and a parabo 1 Ic curve in each dimension fitted, which can 
then be used for a quick and accurate location of the extremum. Such a 
procedure Is completely equivalent to linearization of the fitting 
function by means of a Taylor series and subsequent use of linear least 
squares fit techniques to find corrections to the parameters (see 
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Bovlng:ton 1971) • Standard error analysis techn I qnes appropr la te to 
linear problems may be applied. Such a method vorks well, however, only 
so long* as the starting: point is already Inside the minimum. In a search 
for the minimum, It Is more efficient to simply increment along the 
perhaps circuitous path of steepest descent of X' * Marq[uordt 

algorithm allows a contlnmiro of variation between these two modes, 
depending on the texture of JO space, by the addition of a self-scaling 
constant to the diagonal terms of the normal matrix. The computer code 
first searches for a minimum by moving along the direction of largest 
negative gradient and then, once inside, solves for the exact location of 
the minimum through parabolic Interpolation. Restrictions on the 
allowable values of the parameters are implemented by simply preventing 
the search from proceeding past the 1 1ml t s tops . 

Once 1ns ide the minimum, where the 1 Inear approxlraa t ion to the 
fitting function applies, both estimated errors for the parame ters and 
inter— para me ter correlation coefficients can be simply calculated from 
the inverse of the normal matrix. Except In cases where one or more 
walls of the minimum are formed by numerical limit stops, so to speak, 
the error bars thus obtained have been “reasonable**. In other words, the 
power contained in the noise in both the sky plane and the (u,v) plane is 
about the same . 

The parameters, however , are of ten very highly corre la ted . 
Correlation coefficients anywhere between nearly zero and well above .8 
have been observed. This means, In general, that the estimated errors on 
the parameters do not provide a complete picture of the uncertainties in 
the map* Such high degrees of correlation, It might be added, are not 
unexpected. Since there is no positional Informa t ioii conta Ined in the 
closure phase, the brightness d Is tr Ibut ion can “slide around*' on the 
gr id, to the extent that it will reproduce , and s t i 1 1 reproduce all of 
the data equally well. An even more fundamental problem der Ives from 
having performed the fit in the sky plane, as opposed to the (u,v) -plane. 
Although the grid points are uncorrelated, according to the sampling 
theorem, for the highest resolutions, as far as the majority of lower 
resolution data points are concerned, flux could equally well be 
redistributed between many of the grid points. Within the context of 
sky-plane griddlng, and model fitting in general, there is no way to 
Improve the parameter Independence. 

Since, therefore, the error bar estimates do not mean very mtich, we 
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ore more or 


less left to our own devices to decide liow much faith to put 
In the computer output. Although the numerous trials that can be 
necessary to find a solution may be a bother, they can help here by 
putting a handle of sorts on the uniqueness problem. If essentially the 
same maps emerge from different Initial grids, one may have some faith 
that the distribution Is well determined by the data. If different maps 
emerge, one should bo wary. Experience has shown that. In most cases of 
decent data distribution, the maps are repeatable. We would like, 
however, to Inspire some degree of confidence by mapping otherwise known 
sources. That Is the subject of the next section. 

D. PERFORMANCE ON TEST DATA 

We need a source whose brightness distribution Is known In advance 
to try the algorithm out on. None exist on the ml 1 1 Isecond o f arc scale 
In the sky, so one has been created In the computer. This source (Figure 
I II. 3), known as "XXXXXXXX" for lack of a better name, consists of three 
circular Gaussian components of various sizes and flux densities. Note 
that since the Gauss Ians never truly go to zero outside of any finite 
grid, a certain amount of error is Introduced by truncating the grid. 

This is, however, realistic, because the same thing would happen for any 
actual source with a halo* 

Source "XXXXXXXX” has been ‘’observed" with the same (u,v)-plane 
coverage that was achieved on 3C 04 (declination +42 degrees) during the 
July 1974 session using all 4 stations (Figure V.2). Goheratlon of the 
"data" in this way Is somewhat unrealistic* In actual data, the 
distribution of closure phase points Is affected by the fringe 
amplitudes, there not being any when the fringes are below the detection 
limit on any baseline. There is not, however, much closure phase data 
missing, and s imp 1 Ic 1 ty d Ic ta tes that the test source data be generated 
in this way. The source was also "observed " wi th the sampling obtained 
on 3C 120 (declination +5 degrees) during the same session, to Illustrate 
the cose of poorer (Ti,v)^plano coverage. 

A number of computer runs were performed. First, the system was 
checked on per fee t da ta , the no-no Ise case* Following that, noise on the 
15 percent level was added to both amp 1 1 tudes and c losure phases. Since 
it is possible, or e ven probab le , that the Onsala calibration scale is In 
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Figure I I 1 . 3 . The computer-genera led test source « co ns 1 s t I iig o f three 

circular gauss Inn componen ts with the hnlf“wldlh shown. Component 
fluxes are given in Jy. 'llie angular scale is In milliseconds of 
arc , 
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error, two runs were made with all amplitudes on the Onsala baselines too 
hlgrli by 20 percent, ami too low by 20 percent ^ An attempt was also made 
to simulate the lack of elevation corrections at Oiisala by skewing the 
calibration scale up and down by multiplication with a parabolic function 
of Onsala hour angle. The amount was 20 percent at hour angle 6 hrs . 
These maps of source XXXXXXXX are presented in Figure 111.4. It Is 
fairly clear that errors on this (reasonable, in practice) 1 eve I have 
little effect on the appearance of the map, although some of the details 
are changed. In general, the algorithm seems to do what one would 
expect; for example, if the scale of the highest resolution Is adjusted 
downward, the map Is appropriately broadened.. 

The cases of perfect and no Isy da ta were also tried on the test 
source at the declination of 3C 120 (Figure II 1. 5). As we see In the 
next chapter, only the U, N, and G stations were used for 3G 120, so only 
these were used In the test as well, the larger and more elongated beam 
obscures the structure. The result resembles some of the actual 3C 120 
observa t ions . 

An Infinite number of different runs on realistic and unrealistic 
test sources could be constructed. In view of the success shown hero, 
however, we quit and proceed to rer^l sources. 





Flffure Hi.4. Ilcsults of mopping tests on the test eonrco ns explained 
In the text. nuscltncs used were IIN , NG, HG, IIS, NS declination 
assumed was +42 degrees. 







Chapter IV 


3C 120 


In spite of the poor (u,v) -plane coverage that It generally 
displays for earth i‘otatlon synthesis (being at a declination of only +5 
degrees) » 3C 120 has been a favorite target for numerous VLB! experiments 
(Kellermanu et al. 1971; Shaffer et al# 1972; Kellermann et al* 19<3; 
Shapiro et al. 1973; Wittels et al. 1975; Schlllzzi et al. 1975). 
Outbursts happen fast and often, so that observations bear fruit after a 
small number of years. Like 3C 273 D, 3G 279 (Vhltney et al. 1971), and 
3C 345 (Wittels et al. 1976), It exhibits motion of its amplitude minima, 
and it lends Itself to the same numerous sorts of explanations. It is 
particularly important to distinguish between actual apparent motion and 
a succession of stationary outburs ts tha t mimic such motion. In previous 
experiments, lack of (u,v) -plane coverage and phase data have made it 
Impossible to reliably make the dltlnctlon. The addition of closure 
phase data now opens up that possibility. 


A. 3C 120 IN A NUTSHELL 

3G 120 has a Seyfert-type spectrum (Biirbldge 1967; Sargent 1967; 
Arp 1968), of Class 1 (Kliachlkian and Weedman 1974) with broadened 
permitted lines but narrow forbidden lines. It is remarkable for its 
unusually strong helium lines (Shields, Oke , and Sargent 1972; Shields 
1974) , for its power in th© Infrared (e. g. Rleke and Low 1972) , and for 
its rapid optical (e. g. Klnman 1968) and radio Ce, g* Medd et al. 1972) 
va r 1 a 1 1 o ns . 

Its optical appearance (Arp 1968) is of "a bright s tar 1 Ike niic leus 
wi th spikes of nebulosi ty rad la 1 1 ng f rom the center out Into an 
oval-shaped region". Elec tronographlc observations ( Wa Iker , P ike , and 
McGee 1974) show, 1 n add It ion , an extended nebulosity of elliptical 
o lit line, possibly an under lying e 1 1 ip t lea 1 galaxy. 

The radio source is complex (Clark et al. 1968b). It contains one 
extended component with a normal rad to spec trum, one that becomes 
optically thick near 20 cm, and one that is optically thick in the 
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Tlie flux density of the variable 


centimeter ranjye and Is variable, 
component has been well monitored at several waveleng^ths (e. 

Kellerinann and Pan I Iny-To th 1960; Medd et al, 1972; Dent and Ko Jo laii 
1972). The behavior seems to be dominated by separate and distinct 
outbursts, each of which fits very well to pred Ic t ions made on the basis 
of an ad iaba t ica 1 ly expanding; cloud of relativistic electrons In a 
magnetic field (van der Laan 1966; Pau 1 1 ny-To t h and Kellermann 1966, 

I960). This apparent harmony between theory and observation has earned 
for 3C 120 a reputation, whether It is deserved or not, as the "'source 
that works". Figure IV. I presents the total source flux as a function of 
time during the period covered by the auasar Patrol, as determined from 
the Go Ids tone total power measurements. 

B. THE DATA 

In general, 3G 120 Is rather well resolved by the ftnosar Patrol 
array. This fact, combined with Its relative weakness, conspired not to 
produce fringes on the Onsala baselines. Therefore, only the data from 
the HNG triplet have been used. Quasar Patrol data exist on these three 
baselines for so ven sess Ions be tween August 1973 and January 1975, and 
before that, for six sessions with HG only. Two HG data sets also exist 
for May and August 1975. 

The most striking feature about this data set is the steady motion 
(see Figure IV. 2), from November 1972 through July 1974, of the minima of 
the visibility function inward in the (u,v) plane. Such behavior can 
only be explained as an effective Increase in the extent of the 
brightness distribution. The orientation of these minima indicate 
extension along a pos i t ion angle o f about 65 degrees . 


C. THE MAPS 

Fits to the spa t laity f i 1 tered data were obtained for all seven 
sessions (Figure IV. 3). Details are given in Table IV. 1, and the maps in 
Figure IV, 4. All of these were obtained by 1 ni t la 1 iza t ion to a s ingle 
point. The unknoim NRAO ca 1 1 bra t ion f ac tor In March 1974 was fixed at 
the crossing point betA^een the HG and NG baselines. Any amb Igui ties tha t 
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Table IV. 1* Summary of 3C 120 solutions 


Sess ion 

Half-power beam 
E/W N/S 

Date 

A 


Grid 

E/V N/S 

\/F 

Peak 

Hap Flux 

Total F! 

Au^ 73 

1.86 

7.60 

89 

23 

6 

4 

1.06 

a 

X 

10^ 

10.1 

12. 1 

Oc t 73 

1.86 

7.50 

83 

28 

9 

4 

1.87 

4 

X 

10^ 

10.8 

11.0 

Mar 74 

1.86 

7.50 

75 

22 

8 

4 

0.73 

8 

X 

10^ 

8.0 

9.7 

May 74 

1.86 

7.54 

53 

17 

. 8 

4 

1.36 

4 

X 

10^ 

7.8 

10.4 

Jul 74 

1.86 

7.44 

89 

29 

8 

4 

5. 12 

4 

X 

10^ 

8.5 

10.2 

Oct 74 

1.88 

7.72 

60 

12 

8 

4 

1.22 

4 

X 

10^ 

6.8 

8.4 

Jan 75 

1.86 

7.68 

74 

24 

7 

4 

0.99 

6 

X 

10^ 

8.4 

10.8 




65T Uouy* 












PAS£L!/VB Ha 



Flgrure IV.O.c* Fits to the filtered IIG baseline amp 1 1 1 udes on 3C 120. 
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Flffure IV.4.a. GC 120 maps for the August niid Oc tober 1 973 , and ^Inrcli 
and May 1974 scss tons. The scale Is I n ml 1 1 Iscconds 
Is at the top and east to the left. The contour Interval Is 
Jy/square msec of arc, or 1.1 x 10 K* 
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Figure IV. 4. b. 3C 120 inaps for the July and October 1974, and January 
1975 sessions. The scale Is In iti I 1 1 Iseconds of arc- Nortli Is at 
the top and east to the left. The contour Interval Is 0*05 
Jy/s<tuare msec of arc, or 1.1 x 10 ^ K. 



were discovered showed significantly poorer fits than the ones shown. 

The contour interval is 0.5 Jy per square msec of arc, or 1.1 x 10 K. 

The models account for all but 1,9 jfc 0.2 dy , which Ic taken to be 
coming from a more extended halo component. The fact that this number is 
apparently constant, rather than propor t Iona 1 to the total flux, speaks 
for Its being truly due to undetected power rather than a calibration 
error in the VLBl system. The halo must be larger than about 0.01, but 
perhaps not much larger. There has been some d 1 f f icul ty in fitting the 
ampl i tude po Ints of the very lowest reso lut Ion. Flux for the ha lo may 
come from the decimeter component in the spectrum or from a superposition 
of past outbursts, or these two may be synonymous. 

A few genera 1 comments concerning the qua 1 1 ty of the maps are 
probably in order. The one for October 1973 is probably not worth very 
much, and the data may be a t fault. That experiment was one of the ones 
severely affected by parity errors and there is alivays the possibility 
that the editing is incomplete. Even for the others, however, the lowest 
contours are irregular and variable In shape. This Is certainly not a 
real effect. Extended structures depend almost solely on the single UN 
baseline, which shows a relatively condensed track in the (u,v) plane. 

Tlie double structure, when it shows, however, is determined by the 
positions and heights of the maxima and minima and is no doubt real. 


D. DISCUSSION 

In cases where the source components are easily separable on the 
maps, individual component fluxes have been computed and added to Figure 
IV. 1. The map was simply divided along an arbitrary line about equally 
spaced between the components and the flux summed on each side. Vliatever 
part of the halo power that appears in the map has therefore been 
included. The uncertainty is due mostly to the difficulty in dividing 
the power properly between the components. A range of approximately plus 
or minus 1 Jy is allowed in the division of flux. This number Is taken 
to be the error bar on the component fluxes. 

Now that the reality of the double structure has been established, 
component separations can be easily and accurately estimated from the 
positions of the minima in the (u,v) plane. The ind Ivldua 1 component 
sizes are not Involved, since. If the minima are observed, there Is no 
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need to perform o least squares fit to all the data* As longr as the 
Internal structure of the components remains poorly resolved, the 
computed separation should correspond to the distance between centroids. 
This separation has been computed (Figure IV. 5) for every session where 
nul Is or minima were observed or inferred to lie Just off the end of a 
(u,v) track. In every case where the position angle of separation can be 
determined. It Is near 65 degrees. Ve have therefore constrained It to 
this value for the less obviovis cases as well. The range of permissible 
values Is determined mainly by the distribution of the data. One 
baseline experiments where minima were visible, namely March and May 
1973, January 1974, and May 1975, have also been Included. One baseline 
data also exist for August 1975; It shows only a single poorly resolved 
component. Although the flux ratios are not well determined by these 
observations, they are consistent with the division suggested by Figure 
IV. 1 . For all sessions through March 1974, the first minimum is the one 
observed. In May and July 1974, It Is the second. There may be a 
systematic difference between estimates made from the first and those 
made from the second null. In that the latter may yle Id larger 
separations. This difference might be explained by internal asymmetry In 
the components, for example. If one or both were brightened toward the 
outer edge. In January 1975, the first minimum was Inferred to lie just 
Inward of the NG track. The Indicated rate of apparent expansion is 
between 1.1 and 1.5 msec of arc per year or about 4 times c. 

From the maps, which unf or t una te ly cover only a short time span on 
the scale of the flux changes, and the associated separation plot we 
believe we can make the following four generalizations. 

First, each total power flare is the result of flux variations In a 
single spatial source component. 

Second, source components move with respect to one another In such 
a way that the overall size of the source Increases with time. As we 
shall see shortly. Internal expansion within each component is also 
sometimes observed. 

Third, the overa 1 1 expans Ion has been observed to be more or less 
continuous and linear over a time period of several flares. The 
expansion Is linear to within the error bars through July 1974. The 
separation In January 1975 is slightly less, but by no means is It as 
small as during 1972 and 73. This suggests that the new components are 
appearing on or near the outer edges of the radiating region rather than 
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Figure IV. o. Equivalent point double separation, in milliseconds of arc, 
for 3G 120 as a function of time. The c ire les are separa t io ns at 
2.8“cra wave length taken from Schilizzl et al. (1975). 



recurrent ly in the center . 

Fourth, the some preferred postlon angle of 65 degrees Is 
maintained through the whole time period. 

Ve begin to suspect the possibility that physical outbursts are 
rarer than the total flux curves would suggest, an "outburst" actually 
Involving the delayed production of a number of total flux peaks. In the 
light of this possibility, let us briefly examine earlier data, \O.DI and 
otherwise. Tlie 1972-73 flares were the strongest ever recorded (see Medd 
et al. 1972t Andrew, unpublished) . Prior to that, nothing comparable had 
occurred since 1967-60. Tlie Intervening period was occupied by a series 
of smaller flares. The earliest 3.0-cro VLBI data. In February and 
November 1971 (Shaffer et al. 1971) on Goldstack only, showed structure 
on the 1 to 2 msec of arc scale, presumably Involving one or more of 
these minor events. The two data sets are not on adjacent flares (see 
Selelstnd 1974) and are considered to be too widely separated In time for 
much sense to be made of the variations. April 1972 data (Ke Hermann et 
al. 1973), also on HG alone, showed a point source of about 3.7 Jy. In 
all probability, this Is the first observation of the rising 1972 flare. 
An unresolved component was also observed on O August 1972 (Ke Hermann et 
al. 1973) at 6.3 Jy and on 29, 30 August 1972 (Shapiro et al. 1973) at 
7.1 Jy. This rate of Increase Is believable considering the total power 
changes, for the rise of the 1972 burst. Meanwhile, however, some 
variation In the amplitudes was observed in June (KcHermann et al. 1973; 
Shapiro et a 1 . 1973) , meaning that either the April and August point 

sources are not the same one or some remnant of the old structure flared 
up enough to produce structure. By October (Kellermann et al. 1973) and 
November (Shapiro et al. 1973), the source Is no longer unresolved. From 
thence proceeds the expansion which we have observed on the three 
base 1 lues . 

We note tlint, on the rising side of the flarer the component In 
(lues t Ion was unresolved. Later observations definitely sho w observab le 
component sizes. S Ize es t inia tes (see Ke I lerinann e t a 1 . 1973; Shapiro et 

al. 1973; Selelstnd 1974; Witte Is et al. 1975) are model dependent v 
apparent 9 uper”re la 1 1 vis t Ic ve loc 1 1 Ics are Indicated. Internal expansion 
ijS also necessary to explain the complete invisibility, on our array, of 
older components. 

Unfortunately, interpre tat ion from the amplitudes of only one 
baseline is difficult. if a minimum Is observed, a separation may be 
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computodf provided that there Is enougfh data to specify which iniulitmm it 
is. if the conipoiioiits of the supposed double are dissimilar, the 
c o itipo ne u t ratio Is resolution de pe nde iit, and an esti mn te of it based on 
the depth of the mliilimiin con be misleading. Care must be token not to 
overs ta te the cone lus ions . 

So me data exist at o the r f re <iuenc le s » Sc h 1 1 I zz I et al. (1 97o ) 
discuss three observations In 1973 and 1974. Although, except for the 
first session, they ore based on relatively few mensurenienls , the 
separation estimates agree well with ours (see Figure IV. 5). The same 
components appear to be observed at 2.B and 3.8“cm wave 1 eng tlis . This Is 
to be expected since the frequency difference is not great. 

Shaffer and Schlllzzl (1975) present observations from May 1974 on 
the Green Bank/Owe us Valley baseline at 18 cth. Their resolution was 
comparable to HN at the present frequency (iiiaximum resolution 18.5 
million wavelengths), and the visibility curves show strikingly similar 
variations. The general size and orientation of the compact source must 
therefore be similar at the two frequencies, this time somewhat more than 
a slight difference. It accounts for between 60 and 70 percent of the 
total flux. Even at this relatively long iravelength, more than half the 
power emanates from the compact radio source. 

In brief summary, we have observed the expansion of a series of 
individual o u t burs t /components that comprise 3C 120. The expansion 
appears to be continuous over a number of total power flares. Also, 
available evidence suggests tliat the compact source has the same or 
similar structure over at least a moderate range of frequency. 


- 49 - 


Chapter V 


3C 84 


Tho 3C 04 r«dlo .ouro. 1. .l»U.r In .In. nnd .p.otrn,. lo 30 120. 

It 1. mncli atr.ng«r, however, so that sensitivity Is less ef a piobl 
It 1. slta.t.d at abent +40 degrees deelln.tlon ,e the (n.vl-plane 
.overage 1. better. Ond 1. eh.nge. , In b.lh t.tal tlu. nnd eerr.lated 
(Inn. nueh n.re slowly. In (sol. onr problem 1. Ihls ease 1. net lo make 
sense of the evolution b.s.d on onr brief -.napslK. 1. " . bn. rather lo see 

the changes at all- 


A. 3C 84 IN A NUTSHELL 

3C 84 (otherwise knevn as Perseus A) Is closer than 3C 12© and 
therefore better studied optically. It Is Identified wl th the Seyfert 
galaxy (Seyfert 1943) NGC 1275. the most prominent member of the Perseus 
cluster of galaxies. Assuming a value of H = 5© Km/s Mpe for the Hubble 
constant. It lies 11© Mpe distant (z=©.©18) Unlike other Seyfert 
galaxies. NGC 1275 Is not obviously a spiral. I t has a bright. 

elllptlcal-looklng central part and numerous absorption clouds and 

emission filaments outside (Durbldge and Burbldge 1965, Lynds 197©). 

These filaments lie all around the galaxy, but the maximum extension 
occurs near 0 degree position angle. However, although the cluster 
contains a large proportion of ellipticals, its most prominent member 
does not appear to be one, the stellar component dlspays an A type 
spectrum (Minkowski 1968). The nuclear spectrum, of Seyfert Class 2 
(Khachlklan and Ueedman 1974), has Une widths of several thousand km/s. 
In addition, many of the f 1 laments show sharp lines at a ve loc I ty 300© 
km/s greater than the redshlft velocity, which normally would suggest 
expnlslonof material from the nucleus. To confuse things, however, this 
High velocity system may be seen In 2 1 -cm absorp t Ion aga Ins t the nucleus 

(DeYoung. Hoberts and Sas law 1973). The nuclear optical continuum 
appears to contain a nonthermal component ( Oke 1968), and optical 
polarization has been detected (Walker 1968). The nucleus also appears 
strong in the Infrared (e. g. Rleke and Low 1972) . X-ray emlss Ion from 


- 50 *- 


the Perseusi cluster Is complex In spatial distribution and appears to 
contain a point source component centered on NCC 1275 (Vo Iff et al. 

1974) . 

Tlie associated radio source 3C B4 consists of a number of spectral 
components of a variety of physical sizes (see e. fs- Ryle and Vlndram 
1968). There Is a component with a normal power law spectrum, and an 
angular size of a few minutes of arc. Tliere Is also a peak in the 
spectrum attributed to a component which becomes optically thick near 800 
raiz. VLBl observations (Clarke ot al. 1969; Purcell quoted by 
Paul lny~To th e t al. 1976) show this component to have structure on the 
scale of 0.05 and t be elongated along a position angle near 0 degrees. 
The third component, the one of concern to this VLBI study, becomes 
optically thick near 10 GHz. The angular structure of this presumably 
nuclear source will be discussed shortly. 

The centimeter wavelength component shows flux variations at all 
frequencies where 1 t can be observed. A par t icu lar ly long his tory exls ts 
near 8 GHz, where He esc hen ( 1961) recorded It at 10 Jy In 1959. Since 
that time, It has sho>m a more or less linear Increase (e. g. Dent 1966; 
Medd et al. 1972; Dent and Kbjolan 1972). During the later of the 
experiments discussed here. In early 1975, the flux density at 7.85 GHz 
exceeded 60 Jy. ITils Is a remarkably large, but also rather slow, flux 
variation. Tliere Is some Ind lea t Ion that the Increase may be Irregular 
(e. g. Medd et al. 1972) with a suggestion of superposed outbursts. Tlie 
Increase Is barely dlscernable In the Goldstone total power measurements 
(Figure V. 1) . 

B. THE DATA 

Inspection of the amplitude data reveals that 3C 84 must be about 
the same angular size as 3C 120 at Its largest. Overall, the source Is 
ve r y we 11 r e s o 1 ve d . A t t he 1 o we s t re so 1 u t ion on HN , about 85 percent of 
the total source flux is correlated, while the amp 1 1 tude on the GS 
baseline never rises above 5 Jy. 

The amplitude curves are filled with a whole range of peaks and 
valleys of d 1 f ferent he ights and depths, In apparent confus ion. If the 
amplitude data are projected onto the (u,v) plane, a reassuring feature 
becomes evident, however; what once was chaos falls reasonably well Into 
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a row of parallel ridges (Figure V.2). This suggests that the brightness 
distribution Is elongated along a position angle near -10 degrees, and 
has a prominent separation of about 5 msec of arc. The maxima and minima 
vary widely In height and depth, however, and the s linpl Ic 1 ty seems to 
break down on the Intercontinental baselines. The components of the 
source must therefore have a complexity of their own. 

Because the amplitudes on the UN baseline reach almost 50 Jy (out 
of a total of about 60), wo should expect to account for most. If not 
all, of the flux In the map. Indeed, the spectrum suggests that small 
components should dominate the power at 7.85 GHz. 

Evidence for changes In the brightness distribution Is dlscernable 
In the data, although It Is slight. The heights of the maxima and depths 
of the minima can be seen to change slowly with time, and some changes 
can be seen In the closure phase. These changes are larger than any 
possible error in the data, Inc lud Ing sys tema t Ic errors, but only, 
slightly. We can expect the changes to be difficult to Isolate In the 
maps. 

C. "LOW RESOLUTION" MAPS 

Because we have a longer time series of data for the UNO triplet 
only, and longer still for HG alone, a series of "low resolution" maps 
using only the.se baselines would be of Interest,. Of these, the HG 
baseline Is the longest: 100 million ravelengths. The half width of the 

beam Is 2.06 msec of arc east/west and 2.61 msec north/south. The second 
inaxlmura In the bas Ic amplitude pattern Is just reached by the HG 
baseline. At this resolution, the brightness distribution should be 
relatively simple. 

The fits to the ( f 1 Itered) da ta appear In Figure V.3 and the maps 
In Figure V, 4. All were obtained with a 5 by 10 point grid oriented due 
north/south. The July 1074 solution was obtained after Initialization 
with several different s ingle po ints . As all the amplitudes and closure 
phases change only slightly with time, all the res t of the solutions were 
initialized with the map output by the July 1974 solution. Most of these 
have also, a t one t Ime or another, been gro^vn from a point source. The 
correction for the scaling problem for the NRAO baselines In March 1974 
was adjusted so that the level of the HN amplitudes is the same as It was 
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Table V. 

1. Summary of low 

reso 

lut ion 

3C 84 

so lut 

ions 




Map Flux 

Total Flux 

Ses s ion 

Half-power beam 
E/ff N/S 

Data 

A 

Grid 

E/W N/S 


Peak 


1972 

1.86 

2.38 

90 

0 

5 

10 

2.44 

8 

X 

10 

10 

40.1 

50.0 

May 73 

1,86 

2.50 

74 

0 

5 

10 

9.46 

6 

X 

10 

47.7 

54.2 

Aug 73 

1.86 

2.52 

66 

11 

5 

10 

2.04 

9 

X 

10 '® 

43.4 

54.4 

Oct 73 

1.86 

2.38 

134 

41 

5 

10 

3.38 

9 

X 

io'° 

47.9 

57.0 

flar 74 

1 . 86 

2.38 

183 

40 

5 

10 

5.96 

1 

X 

10 " 

53.3 

57.9 

Kay 74 

1.86 

2.20 

167 

24 

5 

10 

3.50 

9 

X 

lo'r 

51.2 

57.8 

Jul 74 

1.86 

2.38 

230 

49 

5 

10 

3.45 

9 

X 

lo'" 

53.7 

57.9 

Oc t 74 

1.86 

2.38 

105 

29 

5 

10 

6.97 

9 

X 

10'° 

52.3 

61.5 

Jan 75 

1.86 

2.38 

105 

29 

5 

10 

6.97 

9 

X 

io‘° 

52.3 

61.5 





Flj^ure V.3.n* Fits to the f i I tere<i HN bose 1 i lie ninplitiides for the low 
resolution series of maps of OG 84 > 
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Figure V.3.b, Fils to the f I I lered flG base 1 Ine ainpl I ludcs for the 
resolution series of innps of 3C 04. 
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U'ST houk- 

Figure V.3.C. Fils to the filtered IlG base 1 file amp 1 1 tudes f or the low 
resolution series of maps of 3C 04 . 
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V.3.d. Fits to the HGN closure plinses for the low resolution 
series of maps of 3C 04. 



Figure V. 4.a. Low resolution nia 
1973 observing sessions. 
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for the sessions adjacent in time* 

Vlth the exception of August 1973, all curves pass through the data 
points well and In a consistent manner from session to session. The 
August 1973 solution, which Is poorly determined because of the sparse 
data, shows a good, but different fit. Although the map looks similar, 
the solution cannot be considered very reliable. 

Table V, 1 summarizes the solutions. 

Host of the flux appears to come from the southern end of the 
source. This large blob Is what Is responsible for the great power In 
the HN amplitudes. The expected elongation does appear In the map, at 
the expected position angle. In the NG and HG range, the southern 
component Is well enough resolved so that Its effective flux Is 
comparable to the northern componen t ( s ) , and the observed deep nulls are 

produced , 

As expected. It Is difficult to locate the source of the brightness 
clianges. It appears to be related to a flux Increase In the central or 
northern part of the source, or perhaps both. Tlie southern component 
appears to be unchanged; at least, 1 t reaches the same brightness contour 
In all solutions. Tlie appearance Is of a migration of emission to the 
no rth with, tl me • 

One baseline ( HG) data also exist for three sessions In 1972 (May, 
June, and July), and for May 1973. This Is not enough In Itself for a 
map, but a possible solution has been reconstructed In the following way. 
Since the HN amplitudes changed very little with time when observed, the 
August 1973 HN points were added to these data sets, after being scaled 
with the total flux changes, and the grid was Inl t lal Ized wl th the 
October 1973 solution. One solution ,^s done for the combined 1972 data 
and one for May 1973. These are not considered to be unlgue brightness 

distributions, but rather those which require a minimum of modification 
from the be t ter-de termlned c losure phase solutions. The two maps (Figure 
V.5) seem to continue back In time the trend that iras just noted. The 
northern section of the source has been brightening with time. A further 
data set in 1969 at 5-cm wavelength on HN a lone a Iso sho^ra the 
possibility of a .similar s true Uire . The overa 1 1 size and orientation 
must at least be s imllar . Cohen et al. (1971) present an HG amplitude 
curve at 3.8 cm that could we 1 1 be s 1ml lar to our 1974 curve at 3.8 cms. 
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D. HIGH RESOLUTION MAPS 


Clearly, tlie full resolution of the array is necessary to see 
exactly what Is happening. Because the lines of itiaxima and minima do not 
continue out^vard as entirely straight lines to the length of the Onsala 
baselines, we con expect the structure to be more complex at higher 
resolution. 

Certain problems arise with the addition of Onsala. First, the GS* 
baseline sweeps with remarkable speed across the complexities in the 
(u,v) plane. Tlie density of the data points along the track is unequal 
to the task of following them. Tlie data on GS are sparser than on the 
other Onsala baselines because the Go Ids tone station was generally only 
available for one day during each session. The longest baseline, 
therefore, adds relatively few data points to the solution. It does, 
however, add parameters In abundance," because the grid must be more 
finely spaced to accommodate the resolution. It also worsens the problem 
of Interpararaeter correlation^ For this reason, the GS baseline was not 
included In any of the present solutions, nor any of the closure phase 
data sets containing It. 

Solutions were obtained for the October 1973, and for the March, 

May and July 1974, sessions. The grids consisted of 7 by 14 points, for 
a total of 90 fitted parameters. All were initialized with an 
interpolated version of the lo w reso 1 u t ion May 1974 solution. In 
addition, the May 1974 run was also successfully started from a single 
point source. The details are summarized in Table V.2. The fits to the 
data are presented In Figure V.6 and the maps In Figure V.7. This 
contour interval is 2.2 x 10 K. , twice that of the low resolution maps. 
Not enough data ^ms available from the August 1973 session to attempt a 
so iut ion. 

These fits are not, overall, as good os the low reso lut ion ones. 

Most of the deviation appears the closure phase data Involving Onsala, 
where the data are rather sparse and noisy. Some of the problem may come 
from power lying outside the grid. More likely, however, it is due to 
mutually contradictory data points, perhaps caused by calibration errors 
in Sweden. Recall that the scale is uncorrected for variations in system 
temperature, and the 3G 84 data is collected over a very large range of 
elevations. The source never sets in Sweden and, at the end of the 
track, is skimming the northern horizon in lower culmination. Since the 
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Table V.2. Summary of hierH resolution 3C 84 solutions 


Sess ion 

Half— power beam 

Data 

Grid 



E/W 

N/S 

A 


E/V 

M/S 

Oct 73 

1.16 

1.40 

260 

116 

7 

14 

3.36 

Mar 74 

1.16 

1.40 

406 

140 

7 

14 

3.61 

May 74 

1 . 16 

1.40 

352 

102 

7 

14 

2.94 

Jul 74 

1. 16 

1.40 

545 

150 

7 

14 

3.37 


I 

as 


O O 




Peak 


Map Flux 


Total Flux 


1 

X 

10 '• 

48.1 

57.0 

2 

X 

lo" 

55.0 

57.9 

1 

X 

10" 

52.4 

57.8 

2 

X 

lo" 

53.7 

57.9 
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Figure V.6.a. F 1 l.s to the f I 1 lered UN boseline ainpH tudes for llie IiIkH 
rcsolutron 3C 04 maps. 



- 60 - 





BA5£/./A/^ 



</ST‘ Mouf ' 

Figure V.f).!). Fils to the filtered NO bnscllne nniplltudes for the hlgli 
re so lilt Ion 30 04 maps . 
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Fl/?ure V.C.f . PM tfs to tlio IING triplet closure phnses for the hlgrli 
rcRolntlon OC 04 limps. 
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general variations in the high resolution data are reproduced, we can 
expect the main features of the fine structure to bo real. Ve might, 
however, beware of the details. 

These maps will reproduce the trend In the GS amplitudes, but not 
the exact form of the wiggles. Since the wiggles and scatter are almost 
Indistinguishable without better sampling, we can only say that the 
structure Is consistent with GS as well. 

Although the time span Is much shorter for these maps, changes with 
a consistent trend are visible. There appear to be two components to the 
north of the strong one. During this time period, the northernmost 
component brightened while the central one dimmed, creating the 
Impression of power migrating to the north without any actual motion of 
the "hot spots". As before, the brightness of the southern component 
reached nearly the same value in each session. 


E. DISCUSSION 


Conveniently, our perhaps best set of data was spaced not more than 
one week In time ai-my form a four station run at 2.0 cm using the Bonn 
dlsii with Fort Davis, Texas, NRAO, and Owens Valley In the United States 
( Paul Iny-To th e t al, 1976). Tills data set contains no closure phase 
information, but the coverage was relatively good on six baselines out to 
290 million wavelengths. The author’s adaptation of Pau 1 1 iiy~To t h’ s map 
and ours are presented side by side In Figure V.6. The 2.8-cm map has 
been arbitrarily Inverted, to put the power at the proper end of the 
source. It has about twice the resolution that ours does. Like ours, 
PauHiiy-To th’s map is composed of a string of components more or less 
aligned along the same position angle. Beyond that, however, there 
appear to be major differences. It would perhaps be possible to blend 
Paul Iny-To Ih’ s two southern components to form om* one. But our central 
component occupies the space of a large b lank In Pan I lny~ 1 o th’ s map. 

These differences might be caused by rea I s true tiire changes with 
frequency. Al flioiigh the f ro quency d 1 f f erence is not great, the spectral 
region in question is near the synchrotron se 1 f-absorp t Ion peak. If 
different parts of the source became optically thick at different 
frequencies, one might expect rother drastic changes over relatively 
small differences In frequency. 









Alternatively, amblg^ultles caused by the mapping procedures or by 
differences In the (u»v) -plane coverage could be responsible. To test 
tills possibility, our map was used to generate visibility curves on the 
"FOG" baselines. Tlie Bonn baselines are all longer than our maximum 
resolution, so no useful comparison can be made with them. The 
reproduction was good to within a maximum deviation of less than 10 
percent In isolated spots in the (u,v) plane. Tills suggests that our map 
is applicable to the 2.0-cm data as well with only minor changes. There 
is no evidence for large scale changes In s true ture wl th freq[uency In 
this resolution range. This result suggests that the Pauliny-Toth map Is 
not the proper solution and might not fit closure phase data if there 
were any. This cone lus ion is reinforced by the appearance in the VLB 
Array Report (Cohen et al. 1975b), in the course of Its evaluation of 
VLBI modeling problems, of another possible map based on the Paullny-Toth 
data. It looks much like our present results. 

Our mapping procedure has also been confirmed by means of a 
two-dimensional version (Hutton, unpublished) of the Fourier series 
analysis routine developed by Rogers et al. (1974), and by a 
seml-model-f Itt Ing technique developed by Cotton (see Wlttels et al. 

1976). Both were, of course, done on the 3.0-cin data. In addition, a 
map for early 1973, again at 2.8 cm. obtained by a method developed by 
Fort and Yee (1976) is presented in the VLB Array Report. It also 
closely resembles our map. Tl\e Fort and Yee algorithm has also been 
applied to some of the older, and probably inadequately sampled, Guasar 
Patrol data, producing several possible maps (Fort 1976). The shortest 
baseline used here was HG so that the larger scale structure Is totally 
undoterinlned. Kith the addition of a lUtle imagination an even older 
Canadian model (Legg et al. 1973) fit to some 2.0-cin data can also be 
made to resemble the present results. 

We coiicliule that the construction of a opectral Index map of 3C 04 
should wait until phase data becomes available at other f re cpieuc Ics , or 
until trrlds Initialized from an adequate ly samp led closure-phase solution 
can be least squares fitted to the amplitude data where this Is not 
possible. 

Fringe amplitude data have also been taken on 3C 84 a t 2 cm (Klell 
et al. 1975) and at 10 cm (Shaffer and Schl llzz 1 1975) but no models have 
been fitted. The features of the amplitude curves are consistent with a 
similar structure. 
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None of the source components 


We observe no motion of mills, 
appear to move, during the October 1973 to July 1974 time period, by more 
than about 0.2 msec of arc. The corresponding velocity limit Is 0.3 msec 
of arc per year or 0.5 times c. There appear to be only changes In the 
relative power of the components. They change, at least during this time 
period, however. In such a way that power migrates northward. In the 
direction that might mimic apparent expansion If It were to continue long 
enoug^h 


- 80 - 



Cliapter VI 


CONCLUSIONS 

Ke now draw cone lus Ions , on the performance of both the computer 
algorithm and the radio sources. 

A. THE ALGORITHM 

One Important point can be made In favor of the sky-grid fitting 
operation! It was successful. It succeeded in exchanging information 
about the size of the brightness distribution for whatever information 
might be missing about the visibility function. Apart from the size 
restriction, which could be set arbitrarily, it allowed a completely 
general brightness distribution. The maps were constructed with no 
confusing sldelobes. It was possible, also, to Implement reasonable 
physical restrictions on the maps i for example, that the source be 
positive. Furthermore, the procedvire was less sensitive to initial model 
input than is iisual for VLBI mapping or model iiig. 

The generality, however, required huge numbers of fit parameters, 
and hence huge amounts of computer time. The complicated nature of 
spaces with huge numbers of dimensions further slowed convergence by 
introducing complicated correlations between the parameters. ¥e were, 
therefore, thro^m out on our own, so to speak, for error analysis. Our 
evaluation of the procedure rests, then, on artificially generated test 
sources, on comparison with other methods and data sets, and consistency 
of different results on the same source. 

We must conclude that, when the Information is present in the data, 
and this is a question of coverage, the sky-griddlng algorithm is capable 
of extracting it, given computer time. 


B. THE SOURCES 

In two op t lea 1 ly similar galaxies, we f Ind somewhat d Ivergent rad lo 
behavior. In 30 120, whose total power variations are rapid, we also 
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observe rapid changes In the structure. Each of the total power flares 
Is seen to arise from a distinct structural component and the components 
are seen to separate at super re la t i vis tic ve loc 1 1 ies . This separa t Ion 
has continued more or less linearly over a series of separate flares, 
along virtually the some position angle. In 3C 84, on the other hand, 
the Individual components also change In Intensity but do not move 
measurably. The time scale of the variations Is somewhat slower. 

Let us first consider 3C 120, 

In the first, chronologically, Instance of apparent 
super-re lativlstlc mo t ion ( Will tney e t a 1 , 1971), nine poss lb le 

explanations were suggested, all perhaps "almost as dubious as the like 
number of lives of the proverbial felines". The f Irs t sugges t Ion 
(numbered 1) was that the poor (u,v)-plane coverage was not sufficient to 
specify the brightness distribution, and that the double structure was an 
Illusion. However, It has withstood the test of additional data, some of 
it with closure phase, on five sources: 3C 120, 3C 273B, 3C 279 (Will tney 
et al. 1971; Keller maim et al. 1974; Marandlno, unpublished), BL Lac 
(Clark et al. 1973), and 3C 345 (Cohen et al. 1976; Witte Is et al. 1976). 
Now Is the time to re-^examl ne , maps in hand , the o ther e Ight sugges t Ions . 

Explanations (vii), that the Hubble constant is totally wrong, 
(vlli), that our notions of cosmology are totally wong, and ( ix) , that 
the source components cons 1st of tachyonic ma ter la 1 , are too dras t Ic for 
even this many-fold increase In the amount of data to justify. 

Exp 1 ana tion (vi), tha t the quasar reds hi ft Is not a distance indicator, 
seems less viable when app lied to ga 1 axles . Exp 1 a no t ions (ill), the 
so-called "Christmas tree" model where components blink on and off 
randomly to mimic expaus Ion, is el imina ted by Inspec 1 1 on of the 3C 120 
maps (see also Witte Is et al. 1976 for motion observed In 3C 345). We 
would also expect, if (ill) is the correct explanation, to see apparent 
contractions about half the time. None have been observed. The same 
objection applies to explanation ( v) , that multipath propagation 
phenomena are responsible. It would also seem difficult to produce 
expansion con t inuous over several flares. 

That leaves the cat with only two of her lives, The apparent 
ve loc 1 ty o f separa t Ion could represent, a phase ve loc 1 ty, ra ther than a 
group velocity- Alternatively, actual motion of sources might be 
o c c ur ring. I f t he 1 1 ne o f mo tion ma ke s a s ma 1 1 e no ugh angle ( 1 e s s t ha n 
about 30 degrees In the case of 3C 120’s 4c rate) with the line of 
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sight, special relativistic effects can cause rapid apparent motion 
(e. g. Rees 1966} Cavallere, Morrison, and Sartorl 1971; Behr et al. 

1976). Se lets tad (1974) attempted to synthesize the flux history and the 
existing VLDI data on 3C 120 In the light of a model proposed earlier by 
Ozernoy and Sazonov (1969). Explosive events In the nucleus were assumed 
to repeatedly eject pairs of synchrotron emitting plasmolds In opposite 
directions. The plasmolds expand ad laba t lea 1 ly as they travel outward. 

If the velocities are assumed relativistic, the approaching component (A) 
always appears further advanced in Its evolution than the receding 
component (B), because of the light travel time between the two. He 
fitted a series of parameterized outbursts to the total flux history at 
2.0, 3.0, 4.5, and 6.0 era wave lengtlis , and divided them Into pairs. He 
then confirmed the existence of the multiple source components and 
estimated velocities of separation and expansion from published VLBI 
data. He concluded that the VLBI measurements, sparse though they were, 
confirmed the model and, furthermore, that bo th separa t Ion and expansion 
rates were relativistic. His analysis stopped In mld-1973, nearly at the 
peak of the second member of the strong outburst pair. Our multiple 
baseline data began with this pair. 

The logic of the pairs of ejecta in the Ozernoy and Sazonov model 
is based on the prominence of double structure in many radio sources, 
among them the giant doubles, and on the suggestion of duplicity in the 

total flux curves of a few variable sources (Andrew 1973; Seielstad 

1974). Equally as striking as the double outburst phenomenon in the 
appearance of the 3C 120 total flux curve is the rarity of the really 
strong outbursts. Prior to the ones in 1973, the last were the 
“classics" of 1967-68 that contributed so much to the theory of variable 
radio sources (see Ke.llermann and Paullny-Toth 1960) . The smaller pealts 
in between seem by comparison to be on ly a f ter thoughts . Although our 
observations cover less than even one full cycle, they may contain some 
relevant evidence. Only during the s trong f 1 ares Is there evidence that 
the source and the components are cons is tent ly very small. Both become 

larger and more comp lex dur Ing the later flares. A mod I f lea t ion of tills 

simple ejection model mlglit be made to fit the observations If It were 
allowed that true physical outbursts occur less frequently than 
or Ig Ina 1 1 y proposed and that the ejecta do not fade uniformly as 
ad laba t lea 1 ly expand lug clouds as they pass outward. It Is interesting 
that, although these minor flares do show the characteristic behavior 
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with frequency which suggests optlcnl depth chniiges, Allert Olsent and 
Alter (1976) have polarization data that shows a more complex model to be 
required. In particular, they believe the flares never become completely 
optically thin at these frequencies. They also observe continuity of the 
pos 1 1 I Ion angle of polarization over several total power peaks, although 
a possible change occurred In early 1974 at the time of one of the peaks* 

Any model whclh Involves relativistic ejection of any substantial 
amount of matter involves formidable kinetic energy requirements. Also, 
relativistic ejection generally predicts a large power difference between 
approaching and receding components due to the redshift difference. For 
symmetrical ejection, the case considered by Selelstad, it is not 
possible to produce the 4 c apparent expansion without an apparent flux 
ratio greater thau 10 between the two components at their peaks. 

Clearly, the receding component should not be seen at all. 

Alternatively, one of the observed components might be the central site 
of the ejection of the other. In which case there should be some 
observabe difference between the two, or all of the components that we 
observe are approaching by along different angles to the line of sight. 
Ejection of multiple components in a plane which we are viewing edge on, 
for example, would preserve the preferred position angle. Any ejecta 
moving in directions too far removed from the line of sight might not be 
blueshlfted enough to be observable. Note that this model does not put 
the Earth in a preferred position, because if the line of sight did not 
lie near the plane, no superre la t 1 vis t Ic expansion would be observed. 
Alternatively, a model where stationary source components are 
"turned on” by the passage of a moving excitation front might have some 
advantages. The kinetic energy problems disappear and, since the 
radiating material would not be moving, so do the redshift problems. 
Travel time and phase velocity related phenomena are s 1 1 1 1 a va I lab le to 
produce superre la t Ivls t Ic apparen t mo t Ion. 

A different sort of phase velocity model has been put forward by 
Sanders ( 1974). In this mode 1 , the apparent motion depends on 
observation of different parts of a dipolar magnetic field as synchrotron 
radiation from electrons constrained to move along the field lines le 
beamed toward Earth. A pair of light spots which separate linearly with 
time is predicted for a per fee t dipo lar field. One w^ould expect that the 
complexity, presumably in the magnetic field, needed to explain the 
flares would also produce irregularities In the explans ion velocity. On 
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the contrary, we observe a fairly unlforra expansion over several total 
power peaks. 

A1 though it is not s tr Ic t ly necessary, of course , that 3G 04 end 3C 
120 stem from the same physical mechanism, the VLBl observations do not 
appear to exclude the possibility of their doing so. Although their time 
scales of variation and motion are different, they both do exhibit 
distinct structure components which change in flux. 

The apparent complexity of 3G 04 shoivs mostly on Intercontinental 
baselines. On the short baselines, it could pass as a highly unequal 
double. In other words, if 3G 04 were turned sideways, reduced in flux 
by a factor of five, moved to a declination of +5 degrees, and observed 
with the HNG triplet, it might look quite a bit like some of our 3G 120 
observations. Tlie insensitivity of Onsa la would have disguised any flue 
structure that might have existed In 3G 120. If the two sources are in 
fact physically similar, the high resolution observations of 3C 04 should 
shed some light on the character of the apparent motions. Although each 
component is not observed to move, a series of them might be excited In 
succession to mimic expansion. The differences between the sources 
might, then, boil down to time scale and to the relative spatial scale of 
the I nho mo ge ne 1 t les In the material or magnetic field in the source and 
the excitation mechanism. Or, within the context of an e Jec t Ion mode 1 , 
the difference between the sources might be simply the result of a 
difference In the actual velocities or angles involved. 

Both ejection and phase velocity explanations appear to have 
survived the onslaught of data In modified form* It is clear that 4 lu 
both cases, what has been added is a requirement for complexity. The 
simplest forms of the models are not adequate. 

G. LOOKING AHEAD 

Eventually, VLBI exper Iments wl 1 1 no doubt become even grander in 
scale t ha t t he y n e no w , a nd all but t he 1 1 n 1 es t ho les i n t he ( u , v) p I a ue 
will be filled. The proposed "VLB Network" is already In the planning 
stages for such operations (Cohen et al. 1^75b) . Insofar as better 
coverage always helps, maps can always be better determined. 

It is interesting to speculate, however, wha t wi ud f a 1 Is source 
s true tur e work may reap If and when the p 1 anned techni c a 1 re qu ire men t s 
necessary to measure continental drift are Implemented* The first will 
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be, of course, sensitivity, since wider record liig? bardwldths are planned. 
Work could be extended to weaker quasars and perhaps even some of the 
more ordinary Seyfert galaxy nuclei. Wider synthesized bandwldtlis 
provide, purely and simply, more (u,v)-plone coverage and perhaps some 
spec tra 1 Inf orma t Ion. 

Astrometric data of accuracy useful on the scale of a few 
milliseconds of arc may also result. Closure phase based maps could be 
used to make corrections for source structure to the phase data used In 
the as trome tr Ic/geode t Ic solutions, thereby re faring them to a kno>m 
position in the map. Maps from successive experiments could then be 
superposed In absolute coordinates. Notice that the usefulness of this 
scheme implies the necessity of source structure corrections for good 
base 1 Ine so lut ions . 

It Is apparent from this study that exper iments c lose ly spaced In 
time, approximately every three months, are necessary to follow 
unambiguously the source structure changes in 3C 120. On the other hand, 
data over many Individual outbursts, years in other words, are necessary 
to see any pattern. Spectral Information Is also a must* Tlie 
accumulation of the VT^DI data necessary to specify the mechanism of 
compact radio sources will, of necessity, be a long and slow process. 

The Quasar Patrol has shown, however, the value of a systematic 
monitoring program. 

We shall end with a list of questions which might be attacked with 
the help of very long baseline interferometry, and which might shed some 
light on the problem. 

Do outbursts come in pairs? If so, is there any systematic 
difference In the power, duration, Internal expansion rate, polarization, 
or any other observable property, between the approaching and receding 
components? Do any recurrently appear at the same projected position? 

Or do outbursts come in groups of many? Are the apparent rates and 
directions of component separation, and hence the true rate and angle in 
the ejection model, repeatable in the same source? Can the center of 
ejection, if there Is one, be seen as a radio source? Is there any 
physical motion at all? 
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